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1,  Introduction  ; 

The  problem  of  the  diffraction  of  plane  waves  hy  a  dielectric  wedge 
is  of  considerable  physical  and  mathematical  interest.  It  corresponds  to 
the  passage  of  light  thru  a  prism,  or  the  diffraction  by  a  cold  front  above 
a  plane  earth. 

The  diffraction  of  plane  or  spherical  electromagnetic  ■waves  by  a 
perfectly  cwiducting  wedge  was  studied  by  H,  M,  Uacdonald     ,  who  used  the 
method  of  expansion  In  eigenfunctions  in  cylindrical  coordinates.  More 

recently,  the  diffraction  of  pulses  by  perfectly  conducting  wedges  has  been 

(2) 
treated    ,  using  the  method  of  conical  flows. 

Although  the  geometry  of  the  dielectric  wedge  is  the  same  as  that 
of  the  conducting  wedge,  the  fact  that  the  dielectric  "constant"  is  a 
piecewise  constant  function  of  angle  means  that  the  differential  equations 
and  boundary  conditions  of  the  system  are  not  separable  in  polar  coordinates. 
Hence,  methods  involving  separation  in  polar  coordinates,  such  as  the  use  of 
a  transform  theorem  involving  Bessel  functions  integrated  over  their  order   , 
which  appears  superficially  promising  because  the  boundary  conditions  on  the 
wedge  are  independent  of  the  radius,  may  be  expected  to  fail. 

In  view  of  this  circumstance,  an  approximation  method  has  been 
developed.  If  the  dielectric  constant  within  the  wedge  does  not  differ  too 
greatly  from  that  outside  it,  an  expansion  of  the  solution  in  powers  of  the 
percentage  difference  in  dielectric  constant  should  converge  fairly  rapidly, 
and  so  a  few  terms  of  the  series  should  give  a  good  approximation  to  the 
entire  solution.  In  particular,  we  will  restrict  ourselves  to  considering 

(1)  H,  «,  Macdonald,  "ELectromagnetism",  G.  Bell  1934,  pp.  78  ff, 

(2)  J,  B.  Keller  &  A.  Blank,  "Diffraction  &  Reflection  of  Pulses  by 
Wedges  and  Comers",  N,Y,U.  Hath.  Research  Group  Report  No,  EK-21 

(3)  N,  N.  Lebedev,  "Doklady  Akad,  Nauk.  SSSR  58,  1007  (19ii7) 
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the  first  two  terms  of  the  series.  The  zero  order  term  represents  the  Inci- 
dent field,  and  the  first  order  term  represents  the  diffraction  effects. 
Since  for  a  cold  front  the  percentage  difference  in  dielectric  constant 
actually  is  small,  we  will  restrict  ourselves  to  considering  this  applica- 
tion, SLlthouf^h  the  method  used  Is  not  so  restricted.  Hence,  we  place  the 
wedge  over  a  perfectly  conducting  earth. 

In  Section  II,  we  use  Snell»s  laws  to  obtain  a  "geometrical"  solu- 
tion. In  Sections  III-VI,  the  actual  solution  is  found,  whose  plane-wave 
part  agrees  with  the  solution  obtained  in  Section  II.  The  diffracted  part 
of  the  field  Is  expressed  In  the  form  of  an  outgoing  cylindrical  wave  whose 
angular  pattern  Is  determined.  Finally,  the  details  of  the  application  to 
the  cold  front  are  considered  In  Section  VII. 
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2 .  The  Geometrical  Solution 

V?e  shall  first  discuss  the  plane  rave  or  geometrical  part  of  the 
solution,  which  may  be  obtained  by  using  Snell's  laws  and  neglecting  the 
diffracted  waves. 

V?e  must  first  distinguish  between  the  possible  polarizations.  Let 
us  introduce  a  system  of  coordinates  as  in  Figure  1,  with  the  x-axis  horizontal, 
the  y-axis  vertical,  and  the  z-axls  parallel  to  the  interface.  Thus: 

Figure  1 


Now  ve  wish  to  trept  only  the  case  where  the  waves  are  incident 
normally  on  the  wedge.  Thus  the  direction  of  proDOgPtion  o*"  the  incident  wave 
is  normal  to  the  z-axls.      l^   the  Incident  wave  is  a  uniform  olane  v/ave,  so 
that  it  has  no  variations  of  phase  or  amplitude  in  the  z-dlrection,  the  total 
field  will  also  have  no  variation  in  the  z-direction.  Hence,  v:e   reduce  to  a 
two-dimensional  situation. 

We  may  now  further  separate  the  problem  Into  two  cases.  Any  two- 
dimensional  electromagnetic  field  may  be  separated  into  the  sum  of  two  parts 
in  one  of  which  E  Is  entirely  in  the  z-direction,  and  in  the  other  H  Is  entirely 
in  the  z-direction.  These  will  not  couple  to  each  other,  since  the  Maxwell 
equations  break  up  into  two  groups.  Te   may  accordingly  consider  them  separately. 
Since  the  behavior  of  the  two  situations  does  not  differ  significantly,  we 
shall  only  consider  the  case  where  E  is  in  the  z-directlon,  so  E  =  (0,0,E-,). 


In  this  cs£:e,  pince  we  choose  the  grovmd  to  be  perfectly  con- 
ducting, and  F  is  tangentiel  to  the  ground,  F_  nupt  vanish  for  y  =  0. 
Accordingly,  we  describe  a  field  troveling  over  the  ground  not  by  one  but 
by  two  waves  of  eoual  and  opposite  amplitude,  with  eaunl  hori zonal  propa- 
gation constants,  end  eoual  and  opposite  vertical  propagation  constants,  so: 

(2,1)  E2  =  Eq  exp(ik  x  cos  ^Z*)  [exp(lky  sin^)  -  exp(-iky  sin  (^)]  exp(-i6Jt) 

Here  E  is  the  amnlitude,  k  =  (^/z  -  2n/X   ,  and  <^is  the  angle  from  the 
horizontal.  Graphically: 

Figure  2 


y=o 


Here  we  have  chosen  the  angle  ^ to  represent  propagation  to  the 
right.   If  v/2<   ^  <  TT  ,  the  wave  T/ill  propagate  to  the  left.   The  two  cases 
must  be  considered  separately.  We  now  wish  to  have  this  field  configuration 
impinge  upon  the  wedge. 

When  a  wave  strikes  the  wedge,  it  is  partly  reflected  and  partly 

incident  , 

refracted.  Let  us  first  consider /waves  traveling  to  the  right.  The  re- 
flected wave  will  go  off  at  a  direction  making  eoual  ensles  with  the  normal 
to  the  wedge.   The  refracted  wave  will  be  bent  towards  the  normal,  since 
the  wedge  medium  has  a  higher  index  of  refraction  than  the  atmosphere.   If  the 
direction  of  incidence  and  index  of  refraction  are  suitably  relf^ted,  the  wave 
will  he  bent  so  that  it  is  traveling  downward.   It  will  then  strike  the  earth 
and  be  reflected.  Thus  we  hpve  Figure  3^»  3b: 
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Fig.    3a 


Second 
Refracted 
Wave 


Fig.   J>t>. 


Second 
Refracted 
Wave 


The  direction  of  the  refracted  wave  makes  a  smaller  angle  with  the  ground 
than  does  the  normal  to  the  interface,  whence  It  immediately  follows  that 
the  second  refracted  wave  will  not  strike  the  interface,  and  t\6   other 
waves  will  be  generated. 

The  incident  wave  is  present  in  the  entire  region  to  the  left 
of  the  wedge.  This  is  not  true  of  the  reflected  wave.  Since  the  reflected 
wave  must  be  generated  by  waves  striking  the  interface,  it  will  not  anpeer 
in  regions  to  which  such  waves  cannot  penetrate.  Thus,  drew  a  line  parallel 
to  the  direction  of  the  reflected  wave  and  passing  through  the  origin.  The 
reflected  v/ave  will  then  only  exist  in  the  region  between  this  line  and  the 
Interface.  This  Is  because  the  lowest  point  on  the  interface,  that  is,  the 
origin,  gives  rise  to  a  reflected  wave  along  this  dividing  line,  and  no  waves 
can  appear  beyond.  By  identical  reasoning,  the  second  refracted  wave  appears 
only  in  the  region  bounded  by  the  earth  and  by  a  line  through  the  origin  having 
the  direction  of  the  second  refracted  wave. 

Let  us  now  combine  all  these  results  and  apply  to  the  incident 
configuration.  We  then  have  the  following  diagram  for  the  totcl  field  con- 
figuration. We  denote  incident  waves  by  I,  reflected  waves  by  FL,  refracted 
waves  by  RR,  and  second  refracted  waves  by  S.  Subscripts  shov  which  waves  are 
associated.  The  barrier  lines  are  denoted  by  B,(Fpe  Fi^-ure  ^) . 

The  descriptive  behavior  indicated  in  Figure  U  holds  for  all 
angles  9^  between  zero  and'n/2,  although  the  direction  of  the  various  vaves 
and  of  the  barrier  lines  depends  on  the  angle 4^  .  However,  the  behavior  is 
quite  different  when  the  wave  is  incident  from  the  right.  In  thnt  case,  the 
number  of  waves  produced  depends  on  the  angle  of  the  wedge,  ^'e   shall  dis- 
cuss plctorlally  the  behavior  for  a  single  wave  incident  from  the  right  and 
traveling  horizontally.  Though  such  a  wave  cannot  exist  with  our  choice  of 
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polarization,    it  can  exist  when  H  has  only  a   z-component,    the  nltemative 
polarization.        Because  the  pictures  would  become  much  too  cluttered  if 
we  drew  the  two   incident  waves  and   their  behavior,   as   in  Figure   4,   we   shall 
choose  this  simpler  case  of  just  one  wave. 

We  shall  now  draw  pictures  of  the  reflections  and  refractions  for 
several  ranges  of  wedge  angle.      First  suppose  the  angle cX  is   between      TT /2 
end    tt/3.     Then  we  have  the  behavior  indicated  in  Figure  5: 


I- Incident  \^cjve 
M-Reflected  \A^ave 
SRL-Second  ref/ected  ^a\/e 
5RL-  Ba^•^Ier  of  SRL 


^^^re  5  j^^.  Refracted  wave 

SRR -Second  refracted  w<3\/e 
BSRR-Barrierof  SRR 


The  incident  wave  is  reflected  backward  and  downward  at  the  interface,  so  it 
strikes  the  ground  and  a  second  reflected  wave  is  generated.  This  second 
wave  will  not  strike  the  interface  as  long  as  o*  is  greater  than  Tr/3.  At 

^  =  y\/3,    the  second  wave  is  easily  seen  to  be  parallel  to  the  interface. 
The  refracted  wave,  being  bent  away  from  the  normal,  is  bent  downvi-ard  to 
strike  the  ground  and  produce  a  second  refracted  T^ave.  The  second  reflected 
and  second  refracted  waves  both  only  exist  between  the  earth  end  their  re- 
spective barrier  lines. 

Next,  suppose  the  wedge  angle  is  between  tt/3  and  rr/^.,  xhen  the 
second  reflected  wave  will  strike  the  interface  a  second  time.   There,  it 
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generates  two  more  waves,  both  traveling  to  the  right,  one  a  reflected 
wave,  the  other  refracted,  neither  of  which  strikes  the  ground  again.  This 
Is  shown  in  Figure  6. 

Figure  6 


I  -  Incident  Wave 
RL  -  Reflect ed.  Wave 
SEL  -  Second  Reflected  Wave 
TEL  -  Third  Reflected  Wave 
BTRL  -  Barrier  of  THL 
BE  -  Refracted  Wave 
SRR  -  Second  Refracted  Wave* 
TRR  -  Third  Refracted  Wave 
BSRR  -  Barrier  of  SRR 
-  Barrier  of  TRR 


As  can  be  easily  seen,  each  time  the  wedge  angle  ©<■  passes  through^ 
TT /n,  n  =  2, 3, A,  •••  »  nev?  waves  are  generated.   If  n  is  odd,  the  new  waves 
are  produced  by  a  reflected  wave  striking  the  interface  and  undergoing 
partial  reflection  and  refraction.  If  n  is  even,  a  reflected  wave  will 
■strike  the  ground  and  produce  an  additional  reflected  wave.  Thus,  ifcX  lies 
between  'n/n  and  TT/(n  +  1),  there  will  be  n  waves  present  in  the  right  hand 
region,  not  counting  the  incident  wave,  while  in  the  left  hand  region  there 
will  be  [l/2(n  +  3)]  which  is  the  greatest  integer  contained  in  l/2(n+3). 

VTe  thus  see  that  when  the  wedge  angle  is  small,  the  wave  pattern  is  very 
complicated. 
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The  ampldtudes  of  the  waves  can  also  be  determined.  Suppose  the 
change  in  dielectric  constant  is  small,  so  we  write  z  -  (C,  -  €  )/6   where 
subscript  one  denotes  the  left  hand  side,  subscript  two  the  right  hand  side, 
and  z  Is  positive  by  assumption.  Then  we  may  classify  a  wave  by  the  leading 
power  of  z  in  the  power  series  expansion  of  its  amplitude.  For  7  small,  a 
wave  produced  by  e  reflection  at  the  interface  has  an  amplitude  one  order 
higher  than  the  wave  the  produces  it.  Assume  the  Incident  wave  is  of  order 
unity. 

In  Figure  4.,  waves  RL,  and  RL-  are  of  order  z,  while  vraves  RIL  , 
RR2,  S  ,  S  are  of  order  unity.  In  Figure  5,  waves  RL  and  SRL  are  of  order  z, 
waves  RR  and  SRR  of  order  unity.  In  Figure  6,  waves  RL  and  SRL,  and  TF-R  are  of 
order  z,  TRL  of  order  z  ,  and  waves  RR  and  SRR  of  order  unity.  It  is  clear  how 
the  process  continues. 

In  the  next  section  we  shall  consider  how  the  simple  geometrical 
picture  presented  in  this  section  is  modified  by  diffraction  effects.  We 
shall  discuss  how  the  waves  interfere  and  what  is  the  actual  total  field 
amplitude  at  any  point  In  a  later  section. 
3.  Solution  of  the  Field  Eouations 

In  order  to  find  the  total  flelds^including  the  effects  of 
diffraction,  instead  of  just  the  plane-wave  parts,  we  begin  with  Maxvell»s 
eouations: 

I  ■; 

(3.1)  V  X  E  -  ico/4.'K-   0  i 

y  X  ^+  iu>c'f-  0 

Here  the  time  dependence  exp(-lcot)    is  implied.     We  are  interested  "^n  the  case 
where  "e"  has  only  a  z-component,    so  E     =   (0,0,E).     The  first  set  of  Maxv-ell's 
ecuatlons  then  gives: 
(3.?)  H     =    -^—  ^^ 


X       icoyM  8y 

H     -         1        9E 

y       ico/*  a  X 

H     =  0 
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while  the  second  set  yields  the  wave  equation: 

(3.3)        (y^  4-  k^)  E  =  0  k^  =A)3«-£ 

We  are  concerned  with  the  geometry  of  Figure  .1,  and  E  has  no  z  dependence. 
However,  k  is  different  in  the  two  regions.  Calling  the  left  hand  region  1 
and  the  right  hand  region  2,  we  have  k  =  k-,  in  region  1,  k  =  kp  in  region  2, 
where  by  the  physical  conditions  k.^<  k  . 

In  addition  to  the  wave  eauation  3.3,  we  need  a  set  of  boundary 
conditions.   These  are  determined  by  reauiring  that  the  tangential  electric 
field  shall  vanish  at  the  ground,  which  is  taken  to  be  a  perfect  conductor, 
and  that  the  tangential  components  of  electric  and  magnetic  field  shall  be 
continuous  upon  crossing  the  interface.  Let  us  denote  byO(+  an  angle  very 
slightly  greater  than  or  ,  and  by  o<-an  angle  very  slightly  less  than  cK  . 
The  boundary  conditions  are  then  given  by: 
(3.^;)  E  =  0       y  =  0 

jdE(rt4)=  _aiE(  o<-) 
dQ  dO 

The  solution  Is  then  determined  uniquely  by  the  reaulrement  that  it  have  as 

incoming  part  (2.1),  where  it  will  be  convenient  to  take  E  =  1,  and  the 

remaining  part  of  E  shall  represent  outgoing  waves.  V?e  shall  now  develop  a 

method  of  solution  of  these  equations. 

The  wave  equation  (3.3)  may  be  rewritten  in  the  forms: 

(3.5)       ^2    ^2    . 

^-^^^^■'h^^--   (k^  -  ic^)  S(0)E 

S(0)  =  1      0<0<ol 
=  0     o^<  0  <  TT 
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The  form  (3.5)  will  be  convenient  for  treating  reflected  waves,  the  form 

(3.6)  for  transmitted  waves.  In  esoh  form,  the  expression  on  the  right  hand 
side  is  proportional  to  k|  -  k^  =  p.   Hence,  if  this  parameter  p  is  small, 

we  may  treat  the  right  hand  sides  as  perturbation  terms.  V'e  may  therefore 

write; 

(3.7)  E(x,y)  =^p\{x,y) 

o     " 

It  is  clear  that  the  term  n  =  0  represents  the  total  field  when  the  two 
regions  have  the  same  dielectric  constant,  and  thus  it  equals  the  incident 
field.  For  the  higher  order  termn,  we  substitute  the  expansion  (3.7)  into 
the  wave  eauations  (3.5)  and  (3.6)  and  equate  like  powers  of  p  on  both  sides. 
There  results:  i 

^^•^)  (  ^  +^^  +  1^1  )  \i^,y)   =  -  S(Q)  Vi(x,y) 

n=l,2,-.- 

3x   dy     -^ 

Thus,    the   functions  E^^   satisfy  wave  equations  with   ^'orcing  terms.      The 
boundary  conditions   (3.4-)   also  must  be  satisfied  separately  by  each  of  the 
functions  E  .     The  additional  requirement  that  E^  reprenents  outgoing  waves 
at  Infinity  for  n>0  then  determines  F     uniquely. 

To   solve   the   equations    (3.8)   and   (3.9),   we   introduce  the   Green's 
functions:  •.         . 
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G^    (x,x',y,y')   =  \  [H^^k^VCx-x')^  .   (y-/?)   -  i^'^k^  V(x-x' )S   (y.y')')] 


(3.10)  

G^    (x,x',y.y')   ^  [  H^  (k^V(x-x' )%    (y-y'f^)   .  H^^  (k^V(x-x' )2  ^   (y^yV)]     ' 

These  functions  satisfy  the  differential  equation: 

cvx        o)y  ' 

using  an  obvious  shorthand.   Also,  both  G,  and  (^^   vanish  at  y  =  0  and  behave 
like  outgoing  vaves  at  Infinity, 

Ve  now  shall  use  Green's  theorem  to  solve  (3.8).  The  solution 
of  (3.9)  proceeds  in  an  eouivalent  manner.   Green's  theorem  may  be  stated  in 
the  form: 

dx  dy  {e  (x  ,y  )(7  +k^)G  (x,x  ,y,y  )-G-l(x,x  ,y,y  )  (V  ^k^E  (x  ,y  )J 
o  -'-oo  1  ■'■  In 

o>n        Cn 

Here  c  denotes  a  contour  consisting  of  the  x  axis  and  a  large  semicircle  at 

infinity.  Because  of  the  boundary  conditions,  the  integral  along  the  x  axis 

vanishes,  and  the  outgoing  waves  condition  insures  the  vanishing  of  the  integral 

over  the  semicircle.  Using  the  differential  eouations  for  E  and  Gi  gives: 

n     ^ 

0  -  r  r°°dx'dy'{-E^(x',y')i(x-x')/(y-y')+G3^(x,x*,y,y')S(o')E   (x',y*)J 
Jo  J-oo 


(3.13) 


En(^ 


,y)  =  7J  dx'dy'   E^_^  (x',y')  G^(x,x' ,y,y') 
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Haro  we  have  denoted  the  range  of  Integration  by  the  subscript  2,  Indi- 
eating  it  extends  over  region  2,  It  follows  from  (3.13)  that  if  we  know  E  ,,  we 
can  find  EL  by  an  integration.  Since  E  is  known,  we  can  find  the  v'lriouB  E  'a  , 
in  succession.  However,  the  labor  of  finding  E  for  n  >  1  becomes  prodigious. 
But  it  will  develop  that  E.  displays  the  ssuae  geometrical  behavior  as  the  com- 
plete solution,  save  for  certain  amplitude  coefficients,  which  are  correct  to 
first  order.  Hence,  we  shall  restrict  ourselves  to  the  consideration  of  the  term 
El,  which  will  provide  all  the  information  we  need. 


from  (3.13)  satisfies  the  source-free  wave  eouation  in  region  1  only.   Hence  in 
region  1  it  may  be  expressed  in  terms  of  waves  of  the  correct  Dropagntion 
constant.  However,  in  region  2  the  propag'^tion  constants  will  be  incorrect, 
and  we  will  commit  considerable  errors  if  ^"e  investigate  the  f-ir  fields  in 
region  2  on  the  basis  of  (3.13).   It  is  for  this  purpose  th-;t  we  developed  the 
alternative  representation  characterized  by  (3.9).   The  solution  of  (3.9) 
analogous  to  (3.13)  is: 


yx,y)=-ij  -   •  "    "   • 


(3.U)  E^(x,y)   =  -  J^     dx  dy'      E^_^(x   ,y   )      G^{x,x  y,y') 

This  expression  is   suitable   for  analyzing  the   fields  in  region  2. 

An  alternative  method  of  treatment  consists   in  aonlyinf^  Green's 
theorem  directly  to   the  differential  eruetion   (3.5).      This  will   convert   the 
differ'-'ntial  eouation  into  an  integral  ecuation.      Thus,    by  analysis  virtually 
identical  with  that  used  to  derive    (3.13)    from    (3.B),   we  obtain   the   integral 
equation: 


JJ^  dx'dy'      E(x',y')      G^(x,x' ,y,y')   +  E^(x,y) 


(3.15)  E(x,y)   =  p   Jw/^  dx'dy       E(x',y')      G,  (x,x' ,y,y' )   +  E^(x,y) 


15  - 


The  second  term  represents  the  Incident  wave,  that  Is,  thot  due  to  the  sources  at 
infinity.  The  first  term  is  the  analogue  of  a  mass  distribution  in  pot»2ntial 
theory,  or  of  a  distribution  of  simple  sources  in  acoustics.  Such  a  distribution- 
is  knovm  to  give  rise  to  a  field  which  is  continuous  upon  crossing  the  boxmdary 
of  the  region,  and  whose  normal  derivative  is  also  continuous  there.  This  field 
consequently  satisfies  the  conditions  of  the  problem. 

We  may  view  equation  (3.15)  as  constituting  an  integral  equation  for 
the  field  E  when  the  point  x,y  is  in  region  2.  Once  this  equation  is  solved, 
we  may  use  (3.15)  to  evaluate  the  field  at  points  in  region  1,  treating  (3.15)  3S 
the  definition  of  the  field  for  such  points.  To  solve  the  integral  eouation, 
the  method  of  successive  iterations  may  be  employed,  regarding  the  integr&l  terra 
on  the  right  as  a  correction  tsrm.  This  is  equivalent  to  the  Neumann  series 
expansion  of  the  solution.  Thus  we  have  as  a  first  iteration  E  =  E^,  as  a 
second  iteration: 

(3.16)    E(x,y)  =  Eo(x,y)fpJ/dx'dy*Eo(x',y')Gx(x,x',y,y') 

2 

and  so  forth.  This  iteration  scheme,  v^hich  at  each  step  substitutes  the  value 
of  E  from  the  previous  stap,  clearly  leads  to  the  ssme  answer  as  the  per- 
turbation method  described  earlier. 
U.      The  Reflected  Field 

Fe  shall  now  evaluate  the  integral  (3.13)  to  find  the  reflected 
field  to  the  first  approximation.   The  integral  we  need  is: 


u 


,1)  E  =  jJdx'dy'Gi(x,x',y,y')e^V'^°^P'^iy'^^««le-^V'^'"^j 

1     5  *• 
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The  evaluation  of  this  integral  is  based  upon  the  integral  reoresentBtion:^  ' 


U.?) 


^> 


1      ^  '^J-oo  ve^ 


in  which  it  is  assumed  that  k  has  a  positive  imaginary  part.   Bv  means  of 
this  integral,  vre  can  work  out  the  x  ,  y  integrations  of  (U.l) ,   leaving 
the  integration  over  u  to  the  end.  Fe  will  then  obtain  far  field  anproxi- 
mations  by  application  of  the  principle  of  stationary  phase,  "'e  shall  v,rite 


Hence,  we  will  have  to  break  up  (^.1)  into  several  regions.   Let  us  surmose 
the  observation  point  x,y  is  in  the  left  hand  range  with  x  negative.   The 
break-UD  is  indicated  in  Figure  7- 

Figure  J 


Vp  are  dealing  with  the  simolest  possible  case,  since  there  are  only  two  ranges 
of  integration,  those  with  y  ^  y  and  with  y  ^  y. 


(A)      G.  N.  7'atson,  "Theory  of  Bessel  Functions",  Macraillan  19^8,  p.  /J6 

•  The  absolute  value  sign  on  x  is  unnecessary,  but  we  use  the  form  h.?.   to  avoid 

modification  of  the  position  of  certain  subsequently  occurring  sin^jularities 

in  the  coaiplex  u-plane. 
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Thu5  We  have  upon  Fubstltutlng  (A. 2)  into  (ii.l): 


V 


U.3) 


du 


dr^-iky  sinf', 


o    y'cotc< 

■?dy  (^dx'e^^('''-^Y^^''^^'"^^-e'^^^'*"^'f  ^ikx'cos^^riky'sir^-iky'sin^ 
I  y     y'cotW 

The  intef;ration  over  x  end  y  now  involves  only  exDonentlal  functions.  The  algebra 
is  rather  tedious,  so  rather  than  go  thru  the  details  we  shall  only  give  the  result. 
Hence,  lifter  the  integrations  and  simplifications,  we  obtain: 


E, 


1  2 


?  C^-  -iux  ^ 

sino^  I     du  e       I   ivy  j l_ 

^^ij     (ufk  co^)^  /       l^-^k  '^os  {2<(-p) 


u+k  cosi2D(j-p) 


,i(ufk  cos(3P)y  cotP^ 


f  giky  sin( 
/u+k  cos(?o^- 


This  integal  reoresents  a  superposition  of  plane  waves  of  all  directions  of  oro- 
pagation.   Expl'5cit  evaluation  re^^uires  examination  of  the  singularities  of  the 
integrand  in  the  complex  u-olane,  which  isthe  next  order  of  business. 

The  integrand  of  (A. A)  is  vritten  as  the  sum  of  two  terms,  one  in- 
volving exp  ivy  and  the  other  independent  of  v.   The  second  of  these  may  be  evalu- 
ated directly.   This  term  has  a  double  pole  at  u=  -k  cos<pand  simple  poles  et 
u=  -k  cos(2<-f(p)  ,  u=  -k  cos(2t<-<p),  and  no  oth'^-r  singularities  in  the  finite  p^'rt  of 
the  comnlex  u-nlene.   Plso,    because  -xfy  cote(  is  nositive  for  x  negative,  the  path 
of  integration  may  be  closed  bv  an  infinite  semicircle  in  the  unper  half-plane 
and  the  integral  over  the  semicircle  tends  to  zero  as  the  radius  increases.   Hence, 
the  integral  is  eoual  to  the  sum  of  the  residues  at  the  poles  in  the  upper  half- 
plane,  '"e  nor  recall  thf^t  the  imnginnry  p?rt  of  k  is  positive,  ^ov   the  incident 
T«ave  to  travel  to  the  right,  o^f£^,   and  hence  -k  cos<^has  a  negative  imaginary 


IS 


part  and  the  double  pole  is  outside  the  contour  of  integration.  Let  us  now 
suppose  that  the  angle  of  incidence  is  such  that  2d.  -f   and  2  oC  +  f  are  angles 
in  the  second  quadrant.   If  o^  is  fairly  close  to  -  these  conditions  are  satis-  ■ 
fled  for  a  wide  range  of  ^.  The  cosines  of  2  oi.  _  y?  and  2  ^  ■*  f      Tdll  then  be 
negative,  and  the  integrand  will  have  two  poles  in  the  upper  half-plane.   Evaluating 
the  residues  at  these  poles  now  gives: 


^'       ic^ 


(A. 5) 


■^  gik[[cosp-cos   (35['-f)/ cota^-f   sin^yfik   cos(2li'-^)x 

[coaJ?-cos(ai(-f)]  ^ 

1  .ild/cos^cos(3s('f^)|    coIbC-   siitpjy-fik  cos(a:^f^)x 


[cos^-cos(2flf+^)J  2 

Simnlifvinp  the  amolitudes  and  exnonentr  yi'=lds 
(^.6)         E^2=  -  -I 


Ak' 


sin^(c<- 


l eik[  X  cos{2c(-f)  -}- y  sin  {2c(  -<^ 

I ei^[x  cos(?Br'+^  +y  sin  (20^+9)} 


2 
sin  ic{i-q>) 


Upon  consulting:  Figure  \   ve   see  th^.t  this  term  represents  tvo  plane  v.'Hve?,  in 

the  d'rections  of  RL,  and  PUL^  reppectivelv.   T'he  amplitudes  of  there  vaves  thus 

o    2 
have  been  obtained  correct  to  first  order  in  k^  -  k,  .   '"e  thus  hf.ve  the  reflection 

2    2 

coefficients  correct  to  first  order  in  k  -  k 

2    1 


These  may  h<=   dir-^ctly  shovn 
to  be  exactly  the  Fresnel  reflection  coefficients,  apain  to  the  same  order.   Thus, 
the  reflected  nlane  waves  are  the  same  as  those  produced  by  an  infinite  dielectric 
medium  with  the  appropriate  anples  of  incidence.   This  ir  reasonable,  since  the 
nlane  wave  part  of  the  solution  is  mainly  determined  by  the  depree  of  approximation 
of  the  interface  to  a  plane  extending  unv'ard  to  infinity,  rather  th^-n  bv  the  be- 
havior near  the  ground.   We  may  expect  thj^t  the  actual  rpflortion  coefficient  'Aill 
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be  the  same  as  the  Fresnel  coefficient,  correct  to  all  orders.  Horever,  F,^ 
only  renresents  a  mrt  of  the  total  reflected  field.   The  remainder  is  given  by  the 
term  Involving  e  ^  in  (4..4-)»  which  we  shall  call  E   and  will  proceed  to  in- 
vestigate.  E,   has  poles  at  the  same  values  of  u  as  the  poles  of  E.„.   However,  it 
also  has  branch  points  at  u  ■  k  and  u  =  -  k,  because  of  the  appearance  of 


will  remain  a  contribution  from  the  branch  cut  that  mupt  be  put  into  the  u-plane 
to  render  the  integrand  single- valued.  To  evade  this  difficulty,  we  shall  evaluate 
the  integral  approximately  by  use  of  the  principle  of  stationary  phase.  The  re- 
sulting expression  will  reprerent  E-^-^   provided  the  distance  from  the  origin  is 
sufficiently  large,  ind  will  be  fairly  accurate  even  dovm  to  within  a  wave-length 
of  the  origin. 

A  point  of  stationary  phase  exists  where  the  derivative  of  the  exponent 
vanishes.  Writing  x  =  r  cos  0,  y  -  r  sin  0,  where  0  is  an  angle  in  the  second 
quadrant,  the  phase  becomes: 
(4.7)  f(u)  =  r^-u  cos  0  +Vk2-u''  sin  oj 

Different! 8 tiiTg  with  respect  to  u  and  settinr  the  derivative  eoual  to  zero  pives 
the  eouation: 

u  t 

(-C.8)    ,,,       0  =  -  cos  0  -   °  .  sin  0  =  f  (u^) 

A  solution  of  this  is  u  =  -  k  cos  0,  and  it  is  easily  seen  that  this  is  the  only  • 
solution.   Because  0  is  in  the  second  quadrant,  Uq  is  in  the  upper  half-pl^ne.   V'e 
also  have: 

f  (u^)  =  kr 
(i.9) 

f"(u  )  = IL--  

°  k  sin^O 
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V'e  nov  must  consic^er  the  disposition  of  the  poles  and  the  stntionnry 
point  in  the  comolex  u-plane.   First  assvune  7o(\  (p/'^-^>Q,      Then 
-cos  (2<-»^)^  -  cos  (;^-4>)^  -  cos  0,  and  v;e  have  Figure  8: 

Figure  8 


-KCOSO 


In  this  case,  we  may  displace  the  contour  of  integration  from  the  real  axis  to  the 
line  Im  u  =  Im  (-k  cos  0)  without  crossing  any  poles  of  the  integrand.  Calling 
the  factor  multiplying  the  erponentisl g  (u) ,  re  hsve: 


(4.10)  E^ 


^=   /   g(u)e^^(-^u 


0 


Now  f(u)  has  a  strong  maximum  at  u  =  u  .   Hence,  most  of  the  contribution  to  the 
Integral  comes  from  the  neighborhood  of  u  =  u^^,  since  the  exponf'ntiel  oscillates 
rapidly  at  other  points  and  interference  makes  the  total  contribution  avmy   f rom __ 
u  =  Uq  small.  Te   may  accordingly  rrite: 

,co 
U.ll)  Ei^  ^  g  (u^)  eif  K)  I    e^if  (%)(^-%)' 

<.oo  ^" 

where  the  higher  order  terms  have  been  dropped.  This  integral  is  a  standard  form, 
and  we  thus  have: 

U.12)  Eii^g(uQ)  ei^K)  lg^   e"^Z 
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Hence,  in  the  range  2o(  ^.(f>2c(-^>  0,  ve  h^-ve  approximntely 


E  -^  e^^^  e'-^  s J  rTcK sij 


J 


(/;13)^^   VT?  V2^   1,2  (co^cosO)^!   cos(asr-^)-cosO     cos(2^y<p)-cosO 

This  represents  an  outgoing  cylindrical  wave,  vith  a  rnther  complicated 
angulf.r  pattern.   Hence,  for  r  sufficiently  le.rge,  ve  may  omit  this  cylindrical 
wave,  and  the  whole  field  is  eriven  by  (/t.6).   The  ranpie  of  validity  of  this  result 
is  the  region  between  the  barrif-r  line  B-^    (see  Fifnjre  U),  and  the  y-axis,  becauFe 
of  the  assumption  x<o.   However,  it  may  be  shown  by  direct  manipulrtion  that 
(A.^)  represents  E,  throughout  the  range  O^ot*,  which  extends  the  region  of 
validity  up  to  the  wedge  interface.   The  proof  reruirsf  evaluating  the  integral 
(4..1)  with  X  positive,  but  above  the  "edge,  and  differs  only  in  some  algebraic 
details  from  that  used  to  obtain  iU'A) .     Ther-^fore,  we  see  'that  the  exact  theory 
agrees  with  the  geometrical  optics  result. 

The  formula  (-^..1.3)  reDresents  E,,  only  if  we  are  not  too  close  to 
the  barrier  line  0  =  2U-  <Pf    since  there  the  amolitude  of  the  cylindrical  wave  be- 
comes infinite.   This  implies  that  the  field  changes  character  upon  crossing  the 
barrier.  We  shall  investigate  this  point  more  closely  in  a  later  section. 

Next,  Eunnose  that  ?e(  •^(^>  0  >2  0^- <p,  so  the  observation  noint  is  be- 
tween the  barrier  lines.   Then  the  singularities  are  arranged  as  in  Figure  3.. 

Figure  9 
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Not;,    in  order  to  dlsnlace   the  contour  to  pas?   thru   the   5t=tionnry   ooint,   ^e  muct 
cross   the  nole  at  u  =  -k  cos    {20( -<fi) .      Hence,    the  value  of  F^^   is   found   by  adding 
the  residue  at  the  pole   to    U.13) .      This  r^^sidue   is   found   directly   to  be: 


U,u) 


\(P  [C03<f -COsi^-pH     2 


^+  ik  cos(2«V-^)x+ik  sin{^-<j^: 


I e^'^f  ^  cos(?«<-9)fy  sin{7e^-(p)j 


1^     sin^CeC-ip) 

This  residue  will  Just  cancel  the  plane  wave  RL  in  E,/t  Hence,  in  this  region  ve 

have  for  the  total  field  the  plane  wave  RLp  plus  the  cylindrical  vave,  ar;ain  agreeing 

with  geometrical  optics.   This  holds  provided  we  are  not  too  near  P,  and  ^2* 

Finally,  suppose  0>2c<|^.   Tl^en  v^e  see  that  we  must  add  the  residues 

at  both  poles  to  (-^,.13).  The  residue  at  the  pole  u  =  -  k  cos(3^4^)  is: 

_   sin^o<' gikfx  cos(3>rf0  +  y  sin(Ss<+f)] 

(A. 15)  k2[cosf-cos(»^f^)j  ^ 

_  _  JL_    \_ ikfx  cos(2^+<j5»)  +  y  sin(a<+f)] 

^k^     sin2^+,p) 

This  cancels  the  plane  wave  RLo  in  E,o»  so  in  this  range  the  entire  field  is 
given  by  the  cylindrical  wave. 

Consequently,  we  see  that  the  first  order  field  F,  displays  the  beh.-ivlor 
of  the  geometrical  optics  solution,  since  we  have  obtained  a  wave  nicture  for  F-,  which 
is  the  same  as  that  ^iven  by  geometrical  optics.   Moreover,  the  cylindrical  wave 
portion  of  the  total  field,  given  by  (^.13),  shows  thst  the  geometrical  ootics  so- 
lution is  only  valid  at  considerable  distances  from  the  origin,  >^pcaupe  o'l'   the  re- 
ouirements  under  which  the  stationary  phase  analysis  represents  a  pood  aonroximation 
to  the  integral  Enj.  Nevertheless,  the  stationary  phase  method  is  only  off  by  about 
10^  when  we  are  one  wavelength  from  the  origin,  so  the  combinF.tion  of  plane  and 
cylindrical  waves  is  an  accurate  representation  of  the  field  over  most  of  region  1, 
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excent  in  the  Immediate  neiphborhood  of  the  barrier  lines,  where  the  plane  vrave 
is  discontinuous  and  the  cylindrical  wave  is  infinite.   Ve  shall  now  discuss  the 
transmitted  field,  and  then  consider  the  field  near  barriers. 
5.  The  Transmitted  Field 

The  transmitted  field  is  much  more  complicated  to  calculate  than  the 
reflected  field.   Since  the  imaginary  part  of  k  is  taken  to  be  positive,  the  field  E 
is  infinite  at  x  equals  minus  infinity,  and  consequently  the  integrals  do  not 
converge.   In  order  to  obtain  results,  we  must  choose  e  form  for  E^  which  differs 
from  (2.1)  at  infinity,  but  is  close  to  the  correct  incident  field  near  the  inter- 


being  an  incoming  plane  vave  at  infinity,  has  a  change  of  ch'3racter  at  x  =  x  ,  where 

o 

X     is  large  and  negative,    in  that  it  shall  be  an  outgoing  vave  for  x<bc  .     Ye  thus 

choose  for  E^: 

(5,1)  E  =  gik^cos^  I  x-Xq  I  r ^ik^y  sin<p   _^-ikjj   sin<pj 

o 


for  x>x  ,  E  as  chosen  here  differs  from  (?.l)  only  in  a  phase  factor  e~  1   ''^  *, 
and  hence  has  the  same  space,  dependence  as  (?.l).  Ve  now  insert  this  function  into 
(3.1ii)  to  compute  E,  .  Upon  employing  the  integral  representation  (A. 2)  with 
v  =Vk^  -  "  >  we  have  for  E,  : 

E  -  -  1-  r^Jl  dx'dy'  e^^^^/'^'-'^ol+i^^^-^'f^e^^^y-y'^-e^^^^-'^'!)/^^^!^'"^"^-;^^!^  '^) 

(.5.2)    -^       ^ 

The  point  x,y  is  in  region  2,  so  vre  have  to  break  up  the  integral  as  in 

Figure  10, 

*  This  nay  "be  avoided  by  breaking  the  ramge  of  integration  into  that  for  x'  negative 
and  x'  positive,  and  using  to  first  order  tne  known  solution  of  the  problem  of  re- 
flection and  refraction  at  a  plane- interface  to  interpret  the  divergent  x'  negative 
integration. 
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The-re  accordingly  results  for  the  integral  over  x  and  y  in  (5.?),  which  re   shall 
call  I,  upon  omitting  the  lengthy  intervening  algebra,  ' 


c^cos^ 

-ik,  cos(J)XQ+iux+i  (k-j^coa^-u)  co-bX 
^5-^'    +  e  2v  sin^ 


u-k-j^cosp 


-ik,  cosfl>x  4-luxfivy  2 

^  °  2v  sino 


k^sin^-k^-H 


(u-k^cos^)(u-kj^cos(2x'-^)-(k2-k^)sin'W 

-ik]_y  Ein<p 


,2^^,7. 


(u-kj^cop^)  (u-k  cosC^oCf^p)  )-(k^-k'^)sinV 

1 


(u-k-j^cosp)  (u-k-.  cos(2)<'-^))-(k2-kf')£'^noi 


(u-k-j^cos^j)  (u-k^cos(2B<'+^)  )-(k^-k^)  sini^ 
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This  expression  vhen  divided  by  v  represents  the  integrand  of  (5.2), 
from  which  we  are  to  calculate  E,.   Now  suopose  we  set  k  -k_  here;  which  would 
lead  to  no  errors  in  first  order.  The  Integrand  would  then  have  a  double  pole 
at  u=Vj^cos^  .  "  When  complex  integration  methods  are  used  to  evaluate  Ei  in  terms 
of  residues,  the  contribution  from  the  double  pole  would  be  of  order  xe'^^^,  and 
thus  would  represent  a  plane  wave  of  increasing  amplitude.   Such  a  wave,  of  course, 
cannot  he  generated,  whence  it  follows  that  it  is  a  mathematical  fiction  pro- 
duced by  an  inconsistent  set  of  approximations.  We  thus  keep  the  tenn 

2  2    2 
(k2-k2_)sln'£><'  in  the  denominators  of  (5.3),  whence  the  pole  at  u=  k-,cos<^  ap- 
pears to  be  only  of  first  order.   Actually  the  several  terms  of  (5.3)  combine 
in  such  a  way  that  the  integrand  Is  regular  at  u=  k-cos^  .   Accordingly,  we  may 
omit  the  consideration  of  the  pole  u=  k^cos^^  in  the  separate  terms  of  (5.3), 
since  they  will  cancel  upon  combination. 

The  three  terms  of  (5.3)  have  the  following  physical  interpretation. 
The  first  term  cancels  the  incident  wave.   The  second  term  represents  the  refracted 
waves,  find  the  third  term  nrovldes  the  second  refracted  waves  and  also  the  cy- 
lindrical waver  which  reoresent  the  diffracted  part  of  the  field.  Proof  of  these 
statements  will  be  presented  In  the  succeeding  analysis. 

The  two  poles  of  the  first  term  are  located  at  u  =  M  k  -  k^^  sin  (/>\ 
V'e  shall  assume  that  the  imaginary  parts  of  k,  and  kp  are  in  the  same  ratio  as 
the  real  parts,  so  k-j_  and  k2  are  complex  numbers  with  the  seme  phase.   Then  the 
only  pole  of  E,^  in  the  upper  half  plane  is  at  u,^*  Lk^siiKpfVry ,   and  we  may 
close  the  contour  and  evaluate  by  residues,  obtaining: 


h 


(5.0      E^^=  _  £lllf£i?'  e^^^l  '''"''o^  (e^^l^  sin(p_^.i^:^y 


sin^ 


(^H)  2[k^cos^+k2.kJ]' 
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For  small  values  of  p=  k^-kf,    this   Is  approximately: 

^^•^^  ^    ^ ]^     ik^cos^x-x^Y  iky  Eincp     -Ik^y  sin<p) 

^  k2-k2  ^*  "'  ^ 

2     1 

CO 

YTien  this  Is  substituted  into  the  expansion  E=   '\~  p^E  ,   it  just  cancels  the 

o 

term  n=»0."  This  is  proper,  since  the  incident  wave  should  be  removed  to  leave 

only  the  refracted  and  diffracted  fields,  which  are  all  that  are  present  in 
region  2. 

The  integrands  of  K,   end  K   each  have  four  poles.  In  addition, 
E,o  has  a  branch  point  at  u=  kj,  because  of  the  term  involving  v  in  the  exponent. 
Hence,  ve  must  apply  the  stationary  phase  method  to  E,,,  to  obtain  a  good  ap- 
proximation. However,  E^2  ^^Y   ^®  evaluated  explicitly  in  terms  of  residues. 

(5.6)  The  poles  of  E^  are  located  at: 

^21"  k.LCo8c<cos^-<^+sin4<j  kJsin2(o<-^)-fk^-kJ  ]2 
U22=  k^coso^cos^f,p)fsin<5([k^sin2(c^+93)^k^-k^]  ^ 
U23-  k  coso<cosC>(-9?)-sin«<[  k^sin^CoC-^fk^-kM  ^ 
"2a"  kiCop^cos^f^)-sin*c[kjEin2(»c+97)-/k2-k^  j  2 

The  eouare  root  is  so  defined  that  vhen  kp-k  ,  [  k^  sin  ^-^)J   =  k  sin^-^),  and 

similarly  for  [k^sin  6^-M5J)J  ^.  Now  when  k_=  k  ,  we  have  the  limiting  values 

u  =  u  0=  k  cos^,u_-=  k  cos(3s<-<^  ^o/~  k  cos(2s^fa).  Under  the  previous  as- 
2l   22   X        '        ^j        1         •  2^.   1 

siimptions  as  to  the  values  of  the  angles  involved,  u_,  and  yx^  then  have  positive 
imaginary  parts,  while  u^-  and  u_  have  negative  imaginary  p-'.rts.  This  situation 
continues  to  hold  for  k_>k-, ,  since  the  imaginary  p^rts  of  all  the  sauore  roots 
increase.  Hence,  when  we  close  the  contour  with  s  semicircle  in  the  upper  half 
plane,  which  contributes  nothing  because  x-/  cot  JL  is  positive,  (Region  ?.) ,  only 
the  poles  u^^^  aaid  u^^  ^^^^  ^^e  within  the  contour.  The  integral  S  accordincly 
is  the  sum  of  the  residues  at  the  poles  u^  and  u    and  is  thus: 
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(S7) 


F^2  -  \   slrwe   1  o   ^ 


^iu23^>:fii  (k-|_cos9?-U2-]^)  coto<  -^y-^sin^J   y 
,  iu^2^.-fi[{)ii-^coi^-U22)^oiP<    -k,  sin  ^J  y 


2  ,2Ti 


("22  A'^"^^  [^l^^^^f<P)^^2     J 


Ve  shall  r.ov  develop  a  geometric  interpretation  of  the  roots  u^i  and 
.Up^.   Since  U22  differs  from  U21  only  by  changing  the  sipn  of  ^  ,    it  will  be 
sufficient  to  consider  U2-i.  We  shall  show  that  u_  =  k-  cos  B^,  where  &,  is  th 
angle  describing  the  direction  of  propagation  of  the  first  refracted  wave 
(RR  ,  Fig.  U) .   Now  if  p  is  very  small,  the  wave  RR  ,  which  arises  from  an  up- 
going  wove,  (I  ),  will  still  be  upgoing  after  refraction.   As  p  increases,  it 
will  become  horizontal,  and  then  downgoing.   Thus,  the  sign  of  the  angle  B,  wil 
chan/?;e  es  p  Increases.  Fe  shall  define  B  as  positive  for  upgoing  waves  and 


have  Figure  11. 


Fi  gure  1 1 
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By  using  Snell's  lew  and  elementary  geometry,   we  obtain   the   following 
seauence  of  relations: 
(5.8) 


i  = 

^- 

(o^-<?^ 

sin 

i  = 

cos   (o(-^ 

\ 

-  r 

-  (f-<=^) 

cos  B     -  sin  (r+cO 

kj  sin  r  =  k-    sin  i  =  ^]_  cos   (p^-q>) 

kp  cos  B^  =  k2   sin  r  coso(  4  k^  cosTsint^. 

=  k     cos6<     cos(?i-<ff)   f  6ino^[k^  sin^(<<-^i-k^-k^  ] 

=  ^21 
Thus,  the  value  of  u   ,  which  from  (5.7)  is  the  horizontal  component  of  the 

propagation  vector  of  the  wave,  is  the  same  as  that  of  the  wave  RR  ,  determined  by 
geometrical  optics.   A  virtually  identical  proof  shows  thrt  when  the  wave  is  suf- 
ficiently refracted  to  be  bent  downwards,  u_-  still  equals  k^  cos  B^.  7'e  now 
must  discuss  the  coefficient  of  y  in  the  exponent.   Since  the  field  E   satisfies 
the  wave  eauation,  this  coefficient  must  have  the  magnitude  k»\sin  B-^\  .      However, 
the  sign  of  this  quantity  must  be  found.  Fe  have  by  direct  manipulation: 

(5.9)  (k-j^cos^-u^-j^)  coto<'fk^sinp  = 

k  sin<3<cos(D(-<0  -  coso([k^sin^  (o^-^)  f  k^-kj  ]  2         • — 

When  k  =  k  ,  this  expression  eouals  k  sin  <?^.   As  k  inoT'eases,  the  whole  ex- 
pression decreases,  because  of  the  increase  of  the  souare  root,  end  it  reaches  zero 
at  a  value  of  p  equal  to  k  sin(2>^-<^)  sin^/cos  o^  .   For  ls7ger  values  of  p,  (5.9)  be- 
comes negative.   Hence,  the  coefficient  eouals  k  slnB  ,  in  accordance  with  the  con- 
vention that  &|  is  negative  for  down-going  waves.   Thus,  the  field  E,„  represents 
the  refracted  waves.  We  now  expand  the  amolitudes  of  the  waves  in  E   to  first  order 
in  p,  and  obtain  the  simple  result: 
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(5.10) 


2    1 


ikjCx  cos  B,  -^  7  sin  &j^) 


^ikj(x  cos  Bj  +  y  sin  B2) 


Here  the  angle  Bj  Is  that  describing  the  direction  of  propogation  of  the  wave 
RItj,  and  is  negative  for  all  values  of  p. 

Finally,  we  consider  the  field  E,-.  The  integrand  has  poles  in  the 
upper  half  plane  at  U2t  and  U22>  a  branch  point  at  k  ,  which  prevents  closing  the 
contour,  and  a  stationary  phase  point  at  k_  cos  0.  The  stationary  phase  method 
may  be  aoplied  to  evaluate  the  integral,  provided  that  the  contour  may  be  dis- 
placed to  pass  thru  the  stationary  phase  point.  This  requires  investigation  of 
the  disposition  of  the  critical  points  in  the  complex  u-plane.  Now  Im  u  jClm  u^^ 
since  B2  is  clearly  a  larger  (in  magnitude)  angle  than  Rj^.  Hence,  we  have  three 
possibilities.   First,  Im  k2  cos  0  vim  U22»  or  O^Bjl  .  Then  the  contour  may  be 
displaced  to  pass  thru  k  cos  0,  and  evaluation  of  the  integral  yields  the 
cylindrical  wave: 
(5.11) 


=  _  sin  o(   ^^2 


k.sin  0  e-'Vo°°^^   .i(V-p| 


27r      kijcos  0  -  kj^cos  Q>      V    k^r" 


(k^cos  0-u^^)(k2Cos  ©-u^  ) 


(k^cos  ©-U22)(k2Cos  O-Uj^)- 

Hence,  in  the  region  oC>0^/B2|  ,  the  field  is  composed  of  the  refracted 
waves  given  by  (5.10)  and  the  cylindrical  diffracted  wave  given  by  (5.11).  Hext, 
suppose  Im  Upp  <  ln>  k2Cos9<  Im^  U2n->  or  1B21>0>(b  t   Then  we  must  add  the  residue 


ponent  of  propagation  vector  is  lu^  =  k^cosB^,  but  whose  vertical  component  must 


-  30  - 


be  positive,  and  is  thus  kj,|sinB  |.  Upon  expanding  the  amplitude  of  this  wave 

to  first  order  in  p  one  obtains: 

(5.12) 


e-^^l^o^°^f   ik2(x  cos  Bj  f  y  jsin  B^  |  ) 


This  wave  is  thus  the  same  as  the  second  refracted  wave  82-      It  appears  only 
below  the  barrier  line  B2,  where  0  =  |B^3,  exactly  as  given  by  geometrical  optics. 
F\irthermore ,  the  sign  of  the  amplitude  of  this  vrave  is  ODOosite  to  that  of  the 
corresponding  wave  in  (5.10),  thereby  displaying  the  change  in  phase  of  180°  upon 
reflection  from  the  ground. 

Finally,  suppose  Im  U2T_<Iin  k2Cos  0,  or  0<  [&.  \  .  Then,  in  addition  to 
the  residue  at  U22»  we  must  add  the  residue  at  u_  to  the  cylindrical  wave  E  . 
This  residue  is  analogous  to  (5.12),  / 

e'''2 


^^•^^^  _  e-^^l^o^°^^   ik  (x  cos  ^   -f  y/sin  B^D 


This  wave  appears  only  below  the  barrier  line  0  =  fE, | ,  and  corresponds  to  the 
second  refracted  wave  S,  . 

Now  suppose  th^t  p  is  so  small  that  B,  is  oositive,  so  the  wnve  PR 
does  not  strike  the  ground.  Then  we  find  that  for  Q ^  <p ,   the  field  is  given  by 
the  refracted  waves  RJL  and  Pfi^.   For  Q<(p  ,    the  wave  S-,  cancels  the  wave  RR-]_, 
but  the  wave  Sp  is  to  this  approximation  just  equal  to  RP-,,  so  there  is  no  con- 
tradiction.  For  p  sufficiently  large  that  B,  is  negative,  we   have  derived 
exactly  the  field  given  by  geometrical  optics  and  illustrated  in  Figure  4. 

Upon  inserting  the  several  expressions  into  the  expansion  ■  — 
E  =  5:  P""  E„,  and  dropping  the  phase  factor  e'i'^l^o  ^°^^  ''^  ^'^^  ^^^  following 
results: 
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(5.U)  ForcOO^lB^I 

c<  -     iko(^  cos  B,    -^  y  sin  B, )        ik^(x  cos  B^  +  y  sin  Bo) 

h  -  e     d  1         i,_e^       ^         ''tPEc 

For  |B2l>0>|&^l 
E  =  e^^S^^  cos  Bj^  +  y  sin  Bi)_e^^2(^  cos  Bj  -J- y  sin  ^)^q^^2^^  '^°^   ^  +  ^^^^^   ^'^+pE 

For  |B^1>  0 
g  _  gilc2(x  cos  B^  t  y  sin  B^)   ikjCx  cos  B^^  +  y  ^in  B^| ) 
_gik;2(x  cos  Bj  +  y  sin  Bj)   ik2(x  cos  B2  +  y  [sin  B2I  ) 

According  to  (5.L4),  the  transmission  coefficients  of  the  various  waves  are  unity  to 

the  first  order  of  approximation.  Since  the  reflection  coefficients  are  proportional  to 

p,  we  would  expect  transmission  coefficients  behaving  like  1-ap^,  in  agreement  with 

conservation  of  energy.  Thus,  they  should  be  unity  to  first  order  in  p.  Actually 

the  transmission  coefficients  should  be  the  same  as  the  Fresnel  coefficients,  to 

all  orders  of  approximatiDn. 

To  sum  up  our  considerations  on  the  transmitted  field,  we  h<ive  shown  that 

the  field  given  by  the  wave  analysis  agrees  with  that  given  by  geometrical  optics 

in  every   detail,  and  in  addition  T'e  have  found  the  diffracted  field  E  .  All  these 
'  c 

results,  of  course,  are  only  valid  when  we  are  not  too  close  to  the  barrier  lines. 

V'e  shall  novr  develop  the  theory  of  the  behavior  of  the  fields  near  the  barrier 

lines. 

6 .   The  Field  Near  Barrier. Lines 

In  the  neighborhood  of  a  ba-^rier  line,  the  considerations  which  lead  to 
the  cylindrical  vi-aves  have  to  be  modified.  We  recall  that  we  pre  dealing  with  an 
integral  of  the  form: 
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(6.1) 


-^ -no 


where  f  (w)  has  a  zero  at  w  »  w^  and  g(w)  is  suoposed  to  be  slo^rly  varying.   The 

reason  for  the  change  from  u  to  w  will  become  clear  later.   The  major  contribution 

to  the  integral  then  comes  from  the  neighborhood  of  w  =  vi  ,   and  we  replace  g(w)  by 
and  extract  g(w  ) 

g(w  ) /from  under  the  integral  sign.  The  remaining  integral  is  then  evaluated. 

At  the  barrier  lines,  the  stationary  phase  point  of  the  integrand 

coincides  with  a  pole.  We  then  have  g(w  )  =  oo ,  and  the  formula   (A-l?)  certainly 

o 

cannot  be  expected  to  be  valid.  A  modification  of  the  prc^edure,  developed  by 
Pauli   ,  and  Ott   ,  enables  sensible  results  to  be  obtained.  V'e  shall  present 
this  modification  here. 

SuDDOse  the  integrand  of  (6.1)  has  a  pole  at  w  =  w-,  ,  where  w-,  is  close 
to  Wq.  What  we  mean  by  "close"  will  be  made  more  precise  later.  We  than  may 
write  g(w)  =  h(w)/(w-w-j^) ,  where  we  assume  h(w-j^)  to  be  slowly  varying.   This  last 
assumption  is  equivalent  to  stating  that  the  poles  of  g(w)  are  not  close  together, 
so  that  the  maxima  of  the  integrand  do  not  interfere  with  each  other.   We  then 
will  replace  h(w  )   by  h(w  ),  extract  it  from  the  integral,  and  forget  about  it, 
since  it  plays  no  part  in  the  further  calculations.  We  shall  also  assume  thfit 
w  and  w,  are  in  the  upper  half-plane,  with  Im  w-j^>Im  w  .   The  case  Im  r-^tClm  w^ 
will  be  discussed  later.  We  may  then  displace  the  path  to  pass  thru  w  ,  so  f-e  must 
deal  with  an  integral  of  the  form: 


(6.2) 

w-w- 


(5)  V.   Pauli,  Phys.  Rev.  5A,  9^1   (1938) 

(6)  H.  Ott,  Ann.  die  Phys.  ^3,  393  (19^3) 


J-  W-W-.  w 
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The  major  contrihutlon  still  comes  from  the  neighborhood  of  w  .  and 
ve  may  expand  the  exponent,  which  yields  upon  keeping  ouadratic  terms: 
^^•^^  J.  ^if(wjr|lf"(w^)(w-w^)2    _^ 

^  O.  W-Wi 


/_QO  "-(w^-w. 


r-o) 

Under  the  assumption  that  f(v7)  has  a  strong  maximum  at  w  =  w  ,  f  (w  )  is  negative 
and  large  compared  to  unity.   The  integral  may  be  transformed  into: 


=  e^f(-o)[ 
=e^f(-o)/  e-^^   dw 


dw 


-m^ 


(V^'o) 


Here  we  have  introduced  a  dlmensionless  parameter  z,   defined  by  z  =  J_____2__L  f^  ^-^  \. 

V    2     ^  1  o' 
Although  f"('?o)   i?  l«i*KP»  ?•  niay  take  on  any  value  if  W]_-w^  varies.   Hence,  we 

shall  investigate  (6.4.)  for  z  arbitrary.  Ve  note  that  the  imaginary  p'-rt  of  z  is 

positive.   Smell  values  of  z  correspond  to  "closeness"  of  the  Dole  to  the  stationary 

to 

V"e  split  the  integsl  (6.4.)   into  the  two  parts  from  0  to  oo  and  from 

-oo  to  0,  and  write  w  «  -  w  in  the  second.   Then: 

(6.5) 


ooint,  clarifying  the  meaning  of  w-,  close  _. 


Uo      ^-^    Jo  "+2] 

J  o        T^-z2 


)  r  00 

1 


-.,^n-o)   2ze-i-'f°°e-i(^-')  -iw_ 

w^-z^ 
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Fe  introduce  a  real  positive  parameter  t  into  the  exponent,  and  vrite: 

■°°   -it(w2-22)     , 


(6.6) 


F(t,2)  = 

-'  o 


The  integral  in  (6.5)   equals  F(l,z).      Differentiatlnp  F  with  respect  to  t  gives: 

(6.7)  r"^  o     7 

f£=-if      -e^Mw2-z2)       ,, 
©t  Jo 

CO  5 

.    itz2  (         -itw^        , 
=  -  ie    ^'     1         e  dw 

^  o 

The  integral  here  eouals  iVlf  ®   ^  -^ —  »  so  we  have: 

Integrating  with  resoect  to  t  now  gives: 

^  t  2 


F(t,z)  =  F(o,z)  -  I-FK  ^   ^1   ^ 


o   Vr- 

To  evaluate  F(o,z),   we" recall  that  z  has  a   positive   imaginary  part.      The   integral 
(6.6)    for  t  =  o   then  has  only  one  pole   in  the  upper  half-plane,   and   evf^ lusting   the 
residue  at   that  pole  gives   F(o,z)   ='?fi/2z.      The  transformation  T'z^     ~  7  ^   '^^Qn 
throws   (6.9)    into  the  form:  , — 

(^•^Q)  F(t,z)   ^r^  -  Wi   (l-i)   i        J  e'  2  "     dw 

o 

Setting  t  =  1  and  inserting  into  (6.5)  gives: 
(6.U) 


J  .  .''■("o)  TTie-'^'  [  1  -  (1-i)  /"^  '   .'  ^"^  a.] 
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Te   now  vish  to  compare  this  to  the  formula  (^.12)  to  show  the  superiority  of  (6.11). 
It  is  convenient  to  remove  a  factor  27ti    from  both  (^.12)  and  (6.11).   This  is 
because  when  we  consider  the  case  Im  w^  Im  w-,  there  will  be  a  pole-type  contri- 
bution to  the  integral,  and  it  is  desirable  to  make  this  of  unit  amplitude. 

(7) 
The  integral  in  (6.11)  is  in  the  standard  form  of  Freanel' s  integrals 

Te  will  write  u  for  the  ui^per  limit  in  (6.11),  to  apree  with  standard  notation.  Ve 

then  have  the  function: 

^'•"'^  G(u)  =  I  e-'^^i      1  -  (1-i)  f  c(u)  -,  iS(u)]] 

where  C(u)  and  S(u)  are  defined  by: 


S(u)  =  (   sin  -  w2  dw 
J  o 


Upon  expressing  (^..12)  in  terms  of  u,  we  hr,ve  also: 

To  justify  the  notation  G   (u),  we  evaluate  G(u)  for  u  large.  Ve  then  have  the  ap- 
proximations correct  to  order  l/u: 

(6.15)  ■      CCu)-^  H''"^"'  ■  \ 

TT   u 

IT 


S(u) 


cos  2  ^ 


Tf  u 


(7)    E.  Jehnke  znd     F.  Emrie,    "Tables  of  Functions",   Pover  19^5,    n.   34 
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We     then  have  for  G(u): 

(6.16)  G(u)^|  e-'  ?'"'[l  -   (l-i)(^  +  £^  4  |i  -  L^^iLL)] 

•-  ^U  J7-U  -* 

-  i  e-'  I""'  [l  -  [l-l  (1-1)  A.  e'  I'u^]  ] 

Hence,  for  u  large,  G(u)  reduces  to  G  (u) ,  so  (6.11)  and  (A.l?)  agree 

00 

when  u  is  lerge.  Now  large  values  of  u  correspond  to  the  pole  at  w^  not  being 

too  close  to  the  stationary  point  e.t  w  ,  so  we  vould  expect  the  tv'o  fonriules  to 

agree,   ''owever,  (6.11)  is  valid  for  all  values  of  u.   In  particular, 

for  u  =  0,  G(u)  =  \.       Thus,  on  the  barrier  line,  the  field  E  has  the  value 

3  h(w-.)  e    1  .   The  plane  wave  part  of  the  field  F  due  to  the  pole  at  a  has  the 

value  h(w2^)  e  ^''1',  so  the  total  field  on  the  barrier  line  is  .lust  half  the  field 

where  the  olane  v.ave  is  present. 

The  analysis  for  Im  w  >  Im  w.  is  very  similpr.   In  this  case,  v:e  have  to 
add  a  contribution  due  to  the  pole,  and  obtain: 

"•"*         ^r  '^-^  e'^(">a,,-  =  h(.,)a"("l)  +  -1-     (    !ll^  e^f(-'  dw 
2Tt  \\         w-w^  ^   1'  ^  27ri     J     w-w 

-co  ^  —*-       1 

'^o 

•"OO 

The  integral  here  is  evaluated  by  the   same   transfonnptions  that  led   to    (6.11). 
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Calling  u  =  -J^ViffwoV     (w  -w, ) ,    ^6.17)   becomes: 

Vr  ^ 

(6.18)  h(wi)  e^^^"!^  -  h(w^)  e^^^^'o)  G(u) 

For  u  =  o,  this  has  the  value  ^h(w^ )  e    1  ,  and  hence  the  field  is  continuous 
upon  crossing  the  barrier  line.  When  u  is  large,  G(u)  is  small,  and  hence  we 
are  left  with  the  plane  wave  tena  h(w  )  e^^^^l'. 

Graphs  of  the  magnitude  and  phase  of  G(u)  are  given  in  Figures  12  and 
13.  ^e  see  that  the  magnitude  of  G(u)  decreases  steadily  as  u  increases,  and  the 
asymptotic  fom  is  already  Indistinguishable  from  the  exact  formula  when  u  =  2. 
The  phase  tends  rather  slowly  to  ^A,  differing  by  a  ouantlty  of  order  l/u  and 
performing  oscillitions  of  decreasing  magnitude  and  increasing  frequency  as  u 
increases. 

Let  us  now  anply  these  general  considerations  to  a  particular  case.  We 

will  treat  as  a  typical  example  the  reflected  wave  in  the  neighborhood  of  the 

barrier  Bl  (see  Figure  h) ,     We  then  have  w  =  -  k  cos  Q,  w  =  -  k,  cos  (26^  -  <2^  . 

Also  f(w^)  =  k-i^  r,  f(w;,^)  =  k^  r  cosf  0  -(2<p<  -<^)]  ,  h(w^^)  =  sin^c^ /k^(cos  0  -  cosf )' 

h(w  )  »  1/Ak?  sin^{o(-<p),   /f"  (w  )/  =  E— ~  .  We  then  have  for  u: 

^  ^  °     k^  sin^O 

(6.19)  ""^-^IuTq    '^     jcos  0  -  cos  (2<<-<p)| 

Thus,  for  a  fixed  angular  separation  between  0  and  2o(-(p,   u  increases  as  the 
souare  root  of  the  distance  from  the  origin.  Now  the  width  of  the  region  within 
which  G(u)  differs  from  the  asjnmptotic  form  G^u)  is  from  u  =  o  to  u  -  2.  Thus 
for  k^  r  large,  the  angular  width  of  the  transition  region  near  the  barrier  becomes 
small. 
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In  the  transition  region  ve  may  then  approximate  u  by: 

(6.30)  u  =  ^Vv     [Q-(1c(.<p)I 

Thus,  the  total  angular  width  of  the  transition  region  thru  which  E  varies  from 
a  plane  vvave  plus  a  cylindrical  wave  to  Just  a  cylindrical  T.^ave  is  given  by: 

(6.21)  ^0  =4J^=  JL02. 

since  u  must  vary  from  0  to  2  on  both  sides  of  the  barrier. 

This  result  shows  that  for  k  r  large,  the  formula  (4.. 12)  as  derived  by 
stationary  phase  is  valid  everyv.here  except  within  a  range  of  angular  half-width 
^0  =  ^VTf    /V^Tr  around  the  barrier  lines.  Thus,  the  disappearance  of  the  plane 
wave  upon  crossing  the  barrier  takes  place  in  an  angular  range  whose  width  de- 
creases steadily  as  the  distance  from  the  origin  increases.  For  k  r  very  large, 
the  transition  region  is  very  small,  and  a  geometrical  shadow  results.  The  region 
given  by  (6.21)  thus  represents  a  penumbra  region,  betveen  "lit"  and  "unlif 
regions. 

The  curve  u  =  2  is  to  good  approximation  a  parabola  with  the  line 
9  =  2c^-^as  axis.   Only  the  part  of  this  parabola  for  k,r  large  is  of  signifi- 
cance, since  only  then  can  ve  make  the  approximation  that  mopt  of  the  integal 
(6.2)  comes  from  the  neighborhood  of  u  =  u  .   This  limiting  parabola  outside  of 

which  the  plane  and  cylindrical  wave  solution  is  valid  has  appeared  in  the 

(8) 

theory  of  diffraction  by  a  reflecting  half-plane 

To  sum  up,  ve  see  that  the  formula  (6.11)  provides  a  continuous  tran- 
sition for  the  field  E  across  the  barrier  line,  with  a  width  given  by  (6.21). 
On  the  barrier,  the  field  is  reduced  to  just  half  its  "lit  side"  value. 


(8)  B.  Baker  and  E.  Copson,  "Huygen's  Principle",  Oxford  1939,  p.  U5 
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7,  Application  to  Propagation  Through  a  Cold  rront 

The  mathematical  problem  of  the  dielectric  ved^e  is   fairly 
well  realized  in  practice  hy   the  diffraction  of  pl^^ne  waves  by  a  "cold 
front".  This  consists  of  a  mass  of  cold  air  moving  into  a  warmer  region. 
Since  the  cold  air  is  denser,  it  forces  the  warm  air  to  rise  over  it,  and 
the  interface  is  ouite  closely  a  slanting  line.  Actually,  of  course,  the 
air  in  the  immediate  neighborhood  of  the  interface  is  cuite  turbulent. 

To  isolate  the  effects,  we  shall  make  a  number  of  aDproxlmctions. 
First,  we  will  neglect  the.  curvature  of  the  earth,  since  we  are  concerned 
with  fairly  small  regions  of  the  earth's  surface.   Second,  ^ve  shall  assume 
the  earth  to  be  perfectly  conducting,  which  is  well  satisfied  for  most 
freouencies  used  in  practice.  Third,  we  will  neglect  the  turbulent  zone 
near  the  interface.  Under  these  simplifying  assumptions  we  may  use  the 
solution  we  h^ve  obtained  in  Sections  2  -  6  to  discuss  the  propagation  of 
waves  thru  a  cold  front.  The  theory  is  identical  for  a  warm  front,  since 
only  the  direction  of  motion  of  the  air  masses,  which  we  neglect  as  small 
compared  to  the  velocity  of  light,  serves  to  distinguish  between  the  cold 

and  warm  fronts. 

For  an  idea  of  the  order  of  magnitude  of  the  effects  involved, 

we  use  the  well-known  formula'^': 

(7.1)  (n-l)XlO^=  _2.  (p-e-  A8|0e) 

Here  n  is  the  index  of  refraction,  T  the  temperature  In  degrees  Kelvin, 
p  the  pressure  in  millibars,  and  e  the  partial  pressure  of  water  vapor 
in  millibars,  re  also  h'^ve  e  =  .nOl6ps,  where  s  is  practically  eoual  to 


(9)  Burrows  and  Attwood,  "Radio  V^eve  Propagation",  /cademic  Press  19^9 
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the  specific  humidity.      From  these  we  easily  derive: 

C^-^)  dnXlO^  =  -  ^{n-l)xlO^^^{n-l)xlO^'i-^  (-.00l6cs^7.7ds_  ^  V^^^ 

Now  let  us  assume  thnt  one  of  the  two  media  is  a  dry  atmosphere  at 
sea  level  pressure  end  temperature,  so  for  it  p  =  1000  mb,  T  =  300°  K,  and 
s  =  0.  We  shall  take  our  figures  for  the  change  upon  crossing  the  cold  front 
from  a  typical  weather  map,  and  take  dp  =  5  mb,  dT  =  5°  K.  Both  of  thene 
are  somewhat  large.   Also,  "'e  shall  assume  the  atmosphere  is  saturated  on 
the  wet  side,  so  ds  =  1.   This  also  is  an  extreme  figure.   Now  for  these  we 
obtain: 
(7.3)  (n-1)  X  10^  =  263 

dn  -X  10^ '^-.038x263 

T'e  have,  therefore,  a  smell  change  in  dielectric  constant  upon 
crossing  the  cold  front  compared  to  the  dielectric  constant  itself,  and  the 
approximations  of  Sections  2-6  are  valid. 

Te  shall  now  consider  some  practical  aopli cations.   Thus,  consider 
an  airplane  traveling  thru  the  cold  front,  going  from  region  1  to  region  2. 
How  will  its  received  signal  vary  with  position  if  it  flies  at  constant 
height?  V^e  shall  assume  that  the  received  signal  is  proportional  to  the 
magnitude  of  the  field  strength  at  the  airplane,  and  vdll  neglect  variations 
of  the  signal  over  regions  of  space  comparable  to  the  dimensions  of  the 
receiving  antenna,  which  thus  v;ill  be  treated  as  a  point. 

As  long  as  the  plane  is  far  to  the  left,  and  below  the  barrier 
line  B  *  the  received  signal  is  determined  bv  the  incident  wave,  and  thus 
is  proporti onal to  sln(k^y  sina>),  which  is  constant,  thn.t  is,  independent 
•See  Fig.  U. 
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of  X,  at  constant  height  y.  After  crossing  the  barrier  line  Bp,  the  field 

is  the  sum  of  the  incident  field  and  the  reflected  wave  RL  .  Ve  assur.e 

2 

that  k^y  is  large,    so  the  diffracted  cylindrical  wave   is  negligible   compared 

to   the  reflected  plane  wave.      From  eouation   (A. 6)    the    intensity  of  the 

received  signal  is  proportional  to: 

\E|'^2  sin(ky  sin^)4.-E-. _i cosfl^y  sin(2c^+9)-2kX  sin(o^+9)  )Eine<7 

(7.^)  2k  sin{c(-K?)  ' 

Hence,  the  signal  consists  of  a  constant  part  and  a  part  which  oscillates 
with  distance;  that  is,  the  received  signal  is  amplitude  modulated.  Let  the 
velocity  of  the  plane  in  meters/second  be  v.  Then  the  freauency  of  the 
modulation  is: 

(7.5)  fg  =  ^o  •  2  ^  sine<  sin  (<^+<p) 

where  f  is  the  carrier  frecuency  (f  =Ci2/2'7f') .      If  the  carrier  is  in  the 
o  o 

broadcast  band,  f  ^-' 10  cps,  and  the  velocity  is  300  moh,  then  2   f  v/c '^  1  cps 
o  *  o 

and  the  modulation  is  inaudible.  On  the  other  hand,  for  microwave  communi- 
cation, f  /x^x  10°  cps,  we  have  2  f  v/c  "^  3  kc,  so  the  modulation  is  in  the 
audio  region  and  will  be  quite  noticeable  even  though  the  modulation 
coefficient  p/2k^  sin^ (d-i-fp)   is  on  the  order  of  A%.     This  result  is  a  typical 
examole  of  the  inverse  Doppler  effect''   ^  (moving  receiver,  fixed  target),  and 
is  quite  understandable  on  that  basis.   It  also  may  of  course  be  viewed  as  due 
to  interference  between  the  incident  wave  and  the  reflected  wave  RL„. 

When  the  plane  crosses  the  barrier  line  B, ,  the  wave  RL,  is  oresent, 
and  this  produces  further  modulation  of  the  carrier.   In  addition  to  the 
frecuency  f_,  the  received  signal  also  contains  the  freouency  f,,  where:  : 
(7.6) 


:-^  =  f ^  .  2  Z  sinc^  sin  (o<  -(p) 


(10)  E.  J.  Barlow,  Proc.  I.R.E.,^,  3/:^0  (1949) 
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In  case  ^  Is  a  small  angle,  so  the .frenuencies  are  close  together,  the 
modulation  will  display  beats,  where  the  beat  freouency  is  given  by: 

(7.7) 

^B  =  ^o  ^  -  sin  2«i^  sin  (f> 

and  of  course  fo  is  the  difference  betv^een  f,  and  f^. 
a  1  ? 

As  the  plane  continues,  the  same  general  behavior  ensues.  The 

wave  displays  modulation  at  rather  low  frenuencies  between  the  interface  and 

the  barrier  line  B  .  This  modulation  arises  from  interference  between  the 
2 

two  refracted  waves,  which  are  of  the  same  order  of  maenitude.   Thus  modu- 
lation is  maintained  throughout  region  ?.   It  may  be  thought  of  as  arising 
from  the  difference  in  path  length  of  the  directly  received  wave  and  the  wave 
reflected  from  the  grovmd,  which  path  difference  is  prooorticnal  to  the  total 
path  length  and  thus  acts  as  a  modulation.  The  freouency  of  this  modulation 

is: 

2     2 
(7  8)  f  -  ^  f     Z     ^2"*^!  sinO(cos<y  slnC? 

"  ^     o  c       ^2         sini^y  -  <p)  sin(ci -i- <p  ) 

1 

2     2 
Because  of  the  factor  k  -k  ,  the  effective  freouency  may  be  quite  low. 

In  addition  to  producing  modulation,  the  effect  of  refraction  on 

the  transmitted  wave  is  such  as  to  cause  it  to  be  deviated.   Hence,  a  olane 

following  the  beam  to  a  landing  may  not  actually  be  homing.   If  the  plane  is 

flying  toward  the  source,  the  vertical  deviation  of  the  beam  is  in  such  a 

sense  that  the  plane  will  fly  upward  instead  of  horizontally.   Also,  if  the 

beam  is  not  normally  incident  on  the  interface,  there  will  be  a  horizontal 

deviation,  which  will  cause  the  plane  to  follow  a  course  toward  the  apparent 

source  rather  than  the  true  source.   It  should  be  noted  that  all  these 

effects  may  be  masked  by  the  effect  of  the  turbulent  zone  near  the  interface. 


U5 


This  completes  the  discussion  of  the  spatial  behavior  of  the  signal. 
Fe  shall  now  recapitulate  the  results.  Using  MaxT'ell's  eouations  and  suit- 
able approximate  methods  of  calculation,  we  have  obtained  the  fields  yielded 
by  geometrical  optics  and  have  found  the  modifications  on  them  produced  by 
diffraction.  These  results  have  been  expressed  as  combinations  of  waves 
traveling  in  the  proner  directions,  and  the  amplitudes  of  these  waves  have 
been  found  to  first  order  in  the  difference  in  the  propagation  constants  of 
the  media.  Under  the  conditions  of  incidence  which  have  been  treated  here, 
no  new  waves  will  appear  in  the  higher  order  terras  in  the  expansion,  nor  v.ill 
the  directions  of  propagation  be  modified,  although  the  amplitudes  of  course 
will  be  affected.  ^Tien  the  incidence  is  from  region  2,  and  the  angle 
is  small,  new  waves  will  appear  in  higher  approximations,  as  indicated  by 
Figure  6,  but  their  amolitudes  vvill  be  small,  Ve  also  have  discussed  the 
transition  of  the  field  across  the  barrier  lines. 

The  alternative  polarization  (H  horizontal),  may  be  treated  by  the 
same  methods.   The  details  of  the  analysis  are  slightly  more  complicated, 
because  the  boundary  conditions  at  the  interface  are  different,  but  the  results 
should  show  the  same  general  character.   The  same  acnroach  m.ight  also  be 
used  to  treat  such  questions  as  the  wave  theory  for  a  prism  and  the  dif- 
fraction by  bodies  of  polygonal  cross-section. 
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